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Abstract 

We present analytical results for the finite-size scaling in d-dimensional O(N) systems with strong 
anisotropy where the critical exponents (e.g. v\\ and v±) depend on the direction. Prominent examples 
are systems with long-range interactions, decaying with the interparticle distance r as r~ d ~ a with different 
exponents a in corresponding spatial directions, systems with space-"time"a anisotropy near a quantum 
critical point and systems with Lifshitz points. The anisotropic properties involve also the geometry of the 
systems. We consider systems confined to a d-dimensional layer with geometry U n x oo n ; m + n = d and 
periodic boundary conditions across the finite m dimensions. The arising difficulties are avoided using a 
technics of calculations based on the analytical properties of the generalized Mittag-Leffler functions. 

PACS numbers: 05.70.Fh - Phase transitions: general studies; 05.70.Jk - Critical point phenomena. 02.30.Gp - 
Special functions; 
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I. INTRODUCTION 



Anisotropic systems are omnipresent in soft matter and solid state physics. Prominent examples 
are liquid crystals, dipolar-coupled uniaxial ferromagnets, systems with Lifshitz points, systems 
with space -"time" anisotropy near a quantum critical point, systems with long-range interactions 
decaying with the interparticle distance r as r~ d ~ a with different exponents a in corresponding 
spatial directions and some dynamical systems Specific problems arise in the consideration 

of critical phenomena in such systems. In any of the cases the fundamental idea of scaling must to 
be modified in a appropriate way. In terms of the correlation length one can distinguish two type 
of anisotropy: weak and strong. In weakly anisotropic systems, the correlation length has spatially 
dependent amplitude. In the strongly anisotropic systems, in addition the critical exponents (e.g. 
u\\ and u±) depend on the direction. A more general definition based on an anisotropic scale 
covariance of the n-point correlators and different exhaustive examples one can see in ref B2J]. In 
I2I a general approach to scale invariance in infinite volume systems with strong anisotropy has 
been developed. 

The object of the present paper is scaling in finite - size systems. In contrast to the theory 



of finite-size scaling in isotropic systems (see e.g. O, |4|]) and weakly anisotropic systems (see 
e^ Jpland refs. therein) the theory of finite - size scaling in strongly anisotropic systems (see 



e^g J5 land refs . therein) tf 



) is still a field where the lack of results obtained in the framework 
of simplified and analytically tractable models are noticeable. There exist by now quite a few 
examples lEllllKl^l where the predictions of anisotropic finite-size scaling hypothesis have been 
reproduced analytically. In iQ lloi anisotropy appears near a quantum critical point as a result 
of mapping of a "time" dependent problem (in d dimensions) to a "static" problem (in d + 1 
dimensions). In lol it is due to the spatial direction dependence of the interactions . 

Recently J14I a recipe for studying finite-size effects based on some useful properties of the 
generalized Mittag-Leffler functions is suggested. It permits to consider isotropic and some 
strongly anisotropic systems (including long-range quantum systems) on an equal footing. The 
interest in Mittag-Leffler functions has grown up because of their applications in some finite-size 

nnnnn 

scaling problems (see e.g. ||3|,|^ll4j,ll5llla]). The present study (see also 111711 ) is an illustration of 
the rare possibility to handle the final expressions of the scaling equations for strongly anisotropic 
systems analytically. 
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II. THE MODEL 



We restrict our attention to the N-vector spin model defined at the sites of the lattice. The 
Hamiltonian of the model reads 

H = J{x-y)^ y , (2.1) 

where is a classical N-component unit vector defined at site x of the lattice and the spin- 
spin coupling decays with different power-laws in different lattice directions, We assume a d- 
dimensional system with mixed geometry L m x oo™ under periodic boundary conditions in the 
finite dimensions. The interaction between spins enters the expressions of the theory only through 
its Fourier transform. We will consider the following anisotropic small q expansion of the Fourier 
transform of the spin-spin coupling: 

J(q) ~ J(0) + a|||q||| 2CT + ajqj 2 ", (2.2) 

where the first n directions (called "parallel" and denoted by the subscript 1 1) are extended to 
infinity and the remaining m directions (called "transverse" and denoted by _L) are kept finite, with 
m + n = d and a± and an are metric factors and p, a > 0. In finite directions the corresponding 
summations are over the vector qj_ = {q±i, q± m } that takes values in A m defined by q± v = 
2%n u /(aN ) and — (No — l)/2 < n u < (N — l)/2, v — 1, m. In infinite directions the sums 
are substituted with normalized integrals over the corresponding part of the first Brillouin zone 
[— -, ^] n . For our further purposes let us recall that the finite linear dimension L = N a, in the 
continuous limit, means that the lattice spacing a — > and simultaneously N — ► oo. In our 
analysis we accept a± = an = —1/2. Finite-size scaling behavior of such type of systems are 



considered in reference 11311 (with < p, a < 1). In the large-N limit, the theory is solved in 
term of the gap equation for the parameter Ay related with the finite-volume correlation length of 
the system. The bulk system is characterized by a vanishing Aqo, so that the appropriately scaled 
inverse critical temperature 

1 f [7T/a]d dq 
Pc = W? I-*,* |qJ 2 '+|q||l 2ff (2 ' 3) 
is finite whenever the effective dimensionality D = m/p + n/ a is greater than 2. The correspond- 
ing critical exponents v\\ and v± associated with the behavior of the infinite system are 

v± = * ^ - (2.4) 



<t(D-2)' p(D-2)' 

For more details see refJl3l 



III. THE GAP EQUATION FOR THE REFERENCE SYSTEM 



For the system with mixed geometry oo n x L m the gap equation has the form: 

p = — L f 1 "" 1 " y £su (3 1) 

Our analysis will be limited to system below the upper critical dimension D u = 4 and lower 
critical one D\ = 2. 

From the physical point of view, the infinite n-dimensional system, which has a finite size L in 
the remaining m dimensions, can be found in three qualitatively different situations depending on 
the value of -: (i) If 2 < -, then the system is above its lower critical dimension di = 2a and, 
therefore, it exhibits a true critical behavior. A crossover from n-dimensional to D-dimensional 
critical behavior takes place when L —> oo. (ii) In the borderline case of n = 2a, the system is at 
its lower critical dimension and may have only a zero-temperature critical point, (iii) When - < 2, 
the system is below its lower critical dimension and a (D-dimensional) critical behavior appears 
only in the thermodynamic limit L — > oo. 

In the present study we assume that there is no phase transition for finite L, hence n < 2a. For 
n <2a and Ay — > 0, due to the convergence of the integral in (13.11) over d n q\\, one can extend the 
integration over all R n in consistence with the underlying continuum field theory. 

Further the corresponding n - dimensional integral can be presented as 

i s n r° v p^dp 

(2tt)» y |q ± | 2 p + p^ + X v ' 

where S n = 2(7r) n / 2 /r(n/2) is the surface of the n-dimensional unit sphere. With the help of the 
identity: 

v- l dp r(i-;)r(;) i 

rj > a > 0, (3.3) 



t + ^ + |q ± |- rj (t+ |q ± |-)^f ' 

if we choose a = n, r = 2p and rj = 2a, for Eq. d3.2t we end up with the result 



4^ T r^r> 2ff >^ (3-4) 



where 

Sn r(i-£) r(£ 

(2vr)" 2a 



&n,„ = J^h- 1° ^ ■ (3-5) 
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Now the gap equation (13.11) may be presented in the equivalent form 

L qxeA ™ (V + |q±| 2p ) 2CT 

Let us emphasize that one can relate Eq. (l3.6l ) with a fictitious fully finite isotropic m dimensional 
reference system in which the memory of the extended to infinity dimensions and of the anisotropy 
of the system is retained only in the parameter 

Tl 

7 := 1 - — , < 7 < 1, (3.7) 
la 

and in the multiplier A n >cr in front of the sum. 
The normalized m dimensional sum in Eq. (l3.6l) 

wZ£{\v,L):=± E * g. , 2a > n. (3.8) 







can be evaluated with the help of the identity 

1 f°° 

(VT^i ite-»H^EZ m (-X vn (3-9) 

in terms of the generalized Mittag-Leffler function E^^z) (see appendix |A|). If one chooses 
a = p, 7 = 1 — ^ and ?/ = |qj_ | 2 the needed result is 



dxx^El^-Xvx") 



j E ex P(-^ 



geA 1 



7 > 0. (3.10) 



Now let us define 



— o 2 a;) 



aiV, 



iVb/2-l 

aT E ex P 



4tt 2 / 2 x 
' a 2 iV 2 



gGA 1 " l=-N /2 

J I 

and using the approximating formula (5.5) of ref J18ll . we have the expression 



(3.11) 



Qn (x) 



1 



V4 



TlX 



erf 



TlX 



1/2 



2tt 2 x 1 



3 a 



— exp 







[-(;)'* 






y/TTX | 



1 I 00 

) Vexp[-(/aiV ) 2 /4x] 

/7TT ' « 



!=1 



(3.12) 

valid in the large N asymptotic regime. The first and the second terms in the above equation 
are size independent and are precisely the infinite volume limit of Q^ (x). The remainder of the 
calculation involves the insertion of (13.121) into (13.101) . 



IV. FINITE-SIZE SCALING FORM OF THE GAP EQUATION 



In order to illustrate the further calculations we will consider in more details the case m — 1. 
We can represent the right-hand side of Eq. (13.101) as a sum of three terms. 
The first one is given by 



— / dk 

2tt 



i 

H ' w Lix 



dk 



1 



{X v + k 2 P) 1 -£' 



(4.1) 



where the definition of the erf-function 

erf(A-y/x) 

y/Awx 

and the identity (13.91 ) have been used. 
The second term is 

2tt 2 



1 f A 

— / exp(— xk 2 )dk 
2tt y_A 



(4.2) 



3a 



dxx" 1p Ej(-X v x p ) exp 



27P (=) 



7T \2p-l 



3 [Xv + & P Y' +1 



(4.3) 



The third one equals 



dxx^ p - l El lp {-X v x p ) 



TlX 



^exp[-(jiV a) 2 /4x]}>. 



(4.4) 



The first term is exactly the bulk limit W^ 2p (^v, °o). The second one in the considered regime 
Xy — > and a — > is of order O J and can be omitted. 

It is convenient to write the third term, Eq. d4.4ft . in terms of the function (particular case of the 
Jacobi 6 3 function) 

+00 

A(x) = Yl e ~ Xn2 - ( 4 - 5 ) 

I HI— — 00 

and the universal finite-size scaling function Ill9ll 



(2vr; 



2 7 p 



dxx^ p - x El lp 



x h 



(2vr; 



2p 



A m (x) - 1 



7T\ T 

X 



(4.6) 



This can be done with the help of the Poisson transformation formula 

r2~ 



.4 



A{x) = 

and the identity (|A4|>. 

After some algebra the result for the third term is 



(4.7) 



F?, 2p (X v L 2p ) + 



(4.8) 
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Collecting the above results for Eq. (l3.10l) . if m = 1, we obtain 



W? >2p (\ v ,L) = — 

r Z7T 



dk 



2 1P 



(Ay + PP) 1 "^ 



3 [ X y + (I)2P]7 + 1 



7 



— > 0. 
2a 



(4.9) 



In the last term of Eq. (l4. 1 21) apart from the factor L 2 "'^ 1 the intrinsic scaling combination 

y = \ v L 2 ? = (L/U,l? p (4.10) 

emerges, where £x,l is the finite-size transverse correlation length (see Jl^l). The limitation m — 1 
is not principal. If m > 1, in view of Eq. (l3 . 1 2ft . the product [Q No (x)] m in Eq. (l3. 101) contains sums 
of terms of the form 



V4: 



TlX 



erf 



1/2 



exp[- ihaN ) 2 /Ax] 



(4.11) 



8=1 



with 1 < m' < m — 1 and ^ 0, i — 1, m — m'. In such terms the error function erf { ^ x ^ 2 J 
can be replaced by unity, since the exponential function in the right-side of Eq. <!4. Ill) cuts off the 
contribution from values of x 1 ^ 2 <C aN . Note that all the other terms that contain as a multiplier 



(^) 2 x can be estimated. They are of order O (^j^J an d must be omitted in the 



2A1 

3 a eX P 

considered continuum limit. As a result instead of Eq. (l4.9l) we get 

1 f^ m d m k 



L 2 7 p- 



If we introduce the notations 



1 m,2p 



(\vL 2p ) 



1 



(X V L 2 



p 7 



7 = 1 



27 



> 0. 



(4.12) 



Eq. (l3.6l) can be rewritten as 



where 



00 



K-K c oc = Wl 2p (X v , L) - WZ, 2p (0, 00) 



K^(a,p,n,m) = W7 2p (0,oo 



(27T) 



CTk 



(|k| 2p ) 7 



(4.13) 



(4.14) 



(4.15) 
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is the inverse critical temperature (normalized with A„ )0 -) of the "isotropic" bulk system. 
The first term in the right hand side of Eq. (l4. 121) can be presented in the form 

1 f [ ^ ]m . 1 



(Ay + |k| 2 ") 7 



s ^--v r°° x ip -(14- x p v 

K c + _^ X 2 V P 7 / dx-^- 9 / + ) , 7 (4.16) 

2(2vr) m v J x-yp+i-m/2 ^ + xP y 



valid for £j_ )L ^> a. The integral over x converges, provided m > 2^p > m — 2p, and 



x^-'^il + xPf pT(l-f) V 2a 2 P ; - 1 ; 

By substitution Eq. (l4.12l) into Eq. (l4. 14b . taking into account the small-argument expansion 
Eq. (l4.16l) . for the gap equation (13.61) we obtain the scaling form: 



x « -a(n, m; p, a)y D/2 ~ l + F^2 P (y) + i-, 4 > £> > 2, (4.18) 



where 



and 



x = L WP-i)( K - Kc ), y = {L/U,l? p (4.19) 



i sl r< 



m 

n m 



a(n,m;p,a) = - — — ^ — ^— r I 1 - — - — I . (4.20) 

1 P } 2p(27r)-r(l-f) V n 2a 2pJ 

Our model study confirm the phenomenological assumption 111 ill that the finite-size scaling be- 
havior in systems with mixed geometry L m x oo n is governed by the "perpendicular" correlation 
length only. 



V. FINITE-SIZE CORRECTIONS 



Given the gap equation in scaling form, we are now in a position to explore the various finite- 
size corrections. Here, we look at the different regimes: the finite-size scaling regime defined by 
the condition y ~ 1, crossover to the thermodynamic critical behavior y ^> 1, and the regime 
y <C 1. In this section for the sake of simplicity we will consider the important particular case of 
slab geometry, m — 1. 
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A. y ~ 1 



In order to consider the case of y ~ 1 we will use a new representation for F2 2p {y) (see 
appendix IBt 

*wr) - n,m + ^ * * + 2 g ( ^T^iSy . • > ^ 

(5.1) 

and rewrite Eq. (l4. 1 81) in a form suitable for obtaining the shift of the bulk critical temperature. 
Substituting Eq. (l5.1l) in Eq. d4.18l) (m = 1) we obtain the gap equation in the form 



a; 

i 



Therefore, when K — > A'^, simultaneously with L — > oo, in the way prescribed by the equation 

with x = 0(1), the leading-order asymptotic form of Ay is given by 

V - # (5.4) 
where y(a;) is the positive solution of Eq. d5.2l) . Hence, at the critical point x = 0, we obtain 

£±, L = A(n,p,a)L (5.5) 

where A(n, p, a) = l/[y(0)] 1 / 2p is an universal amplitude. In systems with mixed geometry the 
existence of the universal amplitude of the (finite-size) correlation length (i^ on the level of the 



phenomenological scaling has been suggested in 111 ill . Here this qualitative statement is made 
quantitative being a model confirmation of the generalized in Jill Privman-Fisher hypothesis (c.f. 



with Eq.(23) in fllljl) 



When y — > oo, for 2 < D < 4, we may approximate the sum in Eq. d5.2l) by an integral: 

^ (An^r -\y + {±* 2 i 2 YV L_ r , (f ) 2 " p - (i + (f ) 2p Y 7 , = u 2p 

^ (4vr 2 / 2 )^[ ? / + (4vr 2 Z 2 )^ ' (2txz)^p J (£)2 7 p+i/2 (i + (|) 2 p) 7 ' ' U ' 



(5.6) 

With the use of Eqs. d4.17l) and d5.6l) one finds the asymptotic solution of Eq. d5.2l) . which to the 
leading order in x 3> 1 is 

x^-a(n,l;p,a)y D / 2 ^ (5.7) 
i.e. it recovers exactly the familiar bulk hight-temperature result. 



B. y > 1 



The finite-size correction to the bulk critical behavior can be extracted from the asymptotic 
form of the functions Fj a {y) at large argument y ^> 1 (see ful nil ). The result is 



F7,2 P (y) - -y~ J + [2 7 (27r) 2 "C(-2p)] y 

where is Riemann's Zeta function. 

Using Eq. (l5.8l) for the gap equation (13.11) we obtain: 



"(1+7) 



(5.8) 



x 



-a(nA; P ,a)y D / 2 - 1 + [2 7 (27r) 2 ^(-2p)] y~^\ 



(5.9) 



As one can see the finite-size effects governed by the second term in the right hand side of Eq. (l5.8l) 
vary as an algebraic power of the variable y. Since £(— 2p) = for p = k, k - a natural number, 
there are not power low dependent finite size corrections if p — k. The case < p < 1 corresponds 
to the long-range interaction. For p — 1 corresponding to the short-range interaction the result for 
the universal finite-size scaling function is 



-y ' + 



y *e~ 



(5.10) 



[2717(7) J 

that reflects with exponential fall of finite-size corrections in Eq. d5.9l l. And as long as p > 1 and 
if it is not an integer, the power low corrections take place in the case of so-called subleading LR 
interaction but with strong anisotropy. 



C. y«l 

Whenever Eq. (l5.2l) has a solution y 1, use can be made of the asymptotic expansion 

v^, Utt 2 / 2 )^ - \v + (Atv 2 1^pV 

E = - + 2 ^ + °<^ 2 )- (511 ' 



1=1 



In obtaining the first term in the right-hand side of Eq. (l5. 1 ll> the fulfilling of the condition 2^p + 
2p > 1 is used. Taking into account Eq. d5.lll) . in the limit \x — F± 2 (0)\ 1 one finds 

" %-*;U>i'* - <5 ' 12) 

Let us remind, that when the number of infinite dimensions is less than the lower critical dimen- 
sion, the singularities of the bulk thermodynamic functions are rounded and no phase transition 
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occurs in the finite-size system. Nevertheless, one can define a pseudocritical temperature, corre- 
sponding to the position of the smeared singularities of the finite-size thermodynamic functions, 
and study its shift with respect to the bulk value of the critical temperature. In the case under con- 
sideration the first term in the right hand side of Eq. ( 15.21 ) is identified with the shift of the finite-size 
pseudocritical temperature. Actually, for the sake of convenience here we study the quantity K. 
The corresponding result for the pseudocritical K C L is 

Kt - = L~ x Fl 2p (0) , (5.13) 

i.e. the shift critical exponent is A = \/v± in accordance with standard finite-size scaling conjec- 
ture, see Jj]. The coefficient F]_ na (0) can be evaluated analytically as well as numerically for 
different values of the free parameters d, p and 7 using the method developed in reference fl20fl . 



VI. CONCLUSIONS 



The statement, that finite-size scaling in our system with mixed geometry L m x oo n is the 
naturally expected one, is contained in H13I1 where different shape dependent scaling limits have 
been studied. In the present study we are much more interested in the explicit form of the scaling 
equation in different regimes. For this aim an appropriate technics of calculation is developed. We 
show how the mathematical difficulties that arise in the considered anisotropic model with mixed 
geometry can be avoided. First, using the identity (13.31) the problem is effectively reduced to the 
corresponding isotropic one related to a fully finite reference system, Eq. (l3.6l) . A further step is the 
recognition that with the help of the identity (13.91) the appearance of 7 7^ 1 in the summand of the 
gap equation (13.61 ) is not an obstacle for an analytical treatment. Knowledge of the properties of 
the generalized Mittag-Leffler function allows to carry out all calculations analytically. We show 
that though the system is strongly anisotropic, the corresponding gap equation Eq. (l4. 1 81) for the 
intrinsic scaling variable y = AyL 2p has a form similar to the isotropic case with geometry L D (c.f. 
with Eq.(4.115) in We stress that the finite-size L is scaled by the perpendicular correlation 

only. This verified the Privman- Fisher hypothesis for strongly anisotropic systems formulated 

in til- 
Further we conclude that, the finite- size contributions to the thermodynamic behavior decay 
algebraically as a function of L only if < p 7^ k, where k is a natural number. In the case 
p — 1, the finite-size contributions decay exponentially as a function of L. The phenomenon 
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that the so-called subleading terms (in our terminology the term with p > 1) lead to dominant 
finite-size contributions, being unimportant in the bulk limit, was first discussed in ref. ifisl . This 
characteristic feature of the long-range interactions reveals also in our consideration. 
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APPENDIX A: GENERALIZED MITTAG-LEFFLER FUNCTIONS 



I21I (see also 



The generalized Mittag-Leffler functions are defined by the power series 12111 (see also 
) 



fc=0 



(7)* z k 
T(ak + (3) k\ 



a,/?, 7 G C, Re(a) > 0, 



(Al) 



where 



(7)0 = 1, ( 7 ) fc = 7(7 + 1) (7 + 2) • • • (7 + k - 1) = r( * y ) , fe = l,2 ) ... (A2) 

Let us formulate some needed properties of the generalized Mittag-Leffler functions. 
It might be useful to note the relationship 

d 



(A3) 



which follows from the power-series representation. In obtaining Eq. (l4.8l) we have taken into 
account the identity 



dxx^EJ^ (-x p ) = 1, p > 



(A4) 



that follows by integration of Eq. (IA3l) over z from zero to infinity. Next, by subtracting and adding 
l/r(«7) to the function EJ„ /p we obtain 



pry I j.a\ 



y«7 



[1 + Z 



1 + z a Vz a ~i 



(A5) 
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APPENDIX B: DERIVATION OF EQUATION KT\ 



First, we represent the integral in Eq. (l4.6l) as a sum of three terms (m = 1). The first term is 
given by 



Jo 

the second term is 

1 



rW) 



.4 



- 1 



V 



S P,-y(y~' 



(Bl) 



20Fy7-i/2p J Q 
and the third one equals the constant (provided 1 > 2yp) 
1 1 



1 



r(p7) (2tt)^p 



A(x) - 1 



r(P7) 



7T\ 2 



X 



y 



c 

7~l/2p 7>P 



(B2) 



^2p(0). 



(B3) 



Let us now calculate the function S pa (y p ) and the constant C ltP . Making use of the identity dA5l >. 
we represent Eq. dBll) as 

pM?/PJ ^ (4vr 2 / 2 )^[2/ + (47r 2 / 2 )p]7 ' 
To calculate C 7 , p , in Eq. (IB2l) we first write 

1 



(B4) 



t 



-1/2 



7T 



1/2 



dxx 1/2 e tx 



(B5) 



then, by using the identity (IA5I) we take the integral over t, and then 



x p ~< - (1 + x p ) 7 

(1 + £P)7 X P7+V2 : 



(B6) 



i.e. C 7i p 
Eq.CP. 



-am 



1; p, cr). Collecting the results for Eqs. dBll . (IB2I) and (IB3ft for (14.61) we get 
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